Optical properties of the QID multiband models — the transverse equation of motion 

approach 
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The electrodynamic features of the multiband model are examined using the transverse equation 
of motion approach in order to give the explanation of several long-standing problems. It turns 
out that the exact summation of the most singular terms in powers of 1/uj" leads to the total 
optical conductivity which, in the zero-frequency limit, reduces to the results of the Boltzmann 
equation, for both the metallic and semiconducting two-band regime. The detailed calculations 
are carried out for the quasi-one-dimensional (QID) two-band model corresponding to imperfect 
charge-density-wave (CDW) nesting. It is also shown that the present treatment of the impurity 
scattering processes gives the DC conductivity of the ordered CDW state in agreement with the 
experimental observation. Finally, the DC and optical conductivity are calculated numerically for 
a few typical QID cases. 
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I. INTRODUCTION 



Current investigations of the strongly correlated elec- 
tron systems often deal with collective contributions to 
the electrical conductivity (and to the other transport 
coefficients) . This includes in particular the conductivity 
in the charge-density-wave (CDW) or spin-density-wave 
(SDW) ordered states. In most of those studies the elec- 
trical conductivity is determined, theoretically and ex- 
perimentally, relative either to the conductivity in the 
metallic state or with respect to the conductivity in the 
semiconducting state of the pinned CDW/SDW. Obvi- 
ously, the prerequisite to such procedures is the accurate 
evaluation of the conductivity in the metallic state or 
in the state of the pinned CDW/SDW, and this simple, 
yet not completely solved, problem is in the focus of our 
interest here. 

Actually, using the microscopic transverse response 
theory, Lee, Rice and Anderson have found that the 
single-particle optical conductivity in the ordered CDW 
state with the negligible small number of scattering cen- 
ters is given in terms of the semiconducting current- 
current correlation function which describes excitations 
across the gap 0, 0, 0- However, it is also shown 
that for the typical value of the CDW gap and of the 
zero-frequency damping energy (arising from the impu- 
rity scattering processes) their result matches up neither 
the result of the longitudinal response theory 4] nor the 
experimental observation 0, [||. Yet, the longitudinal 
and transverse response have to coincide for fast enough 
(quasi)homogeneous longitudinal fields, as implicit for 
example in the Maxwell equations of the medium, which 
employ only one dielectric function. 

In his textbook pfl, Mahan has further shown that 
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an alternative, so-called force-force correlation function 
method gives a good description of the (high-frequency) 
optical processes in both metallic and semiconducting 
systems, including various excitations within the conduc- 
tion band and across the (pseudo)gaps, but also fails to 
reach correctly the u> — > limit, requiring a specific lo = 
field-theory approach 0, S 13 • 

Also important is the observation that most of the 
transport-coefficient analyses are based on the Boltz- 
mann equations applied to the nearly-free electron mod- 
els jjj llClf , completely neglecting band periodicity in the 
reciprocal space. 

Most of these issues can be settled down using the lon- 
gitudinal response theory with a particular care devoted 
to the continuity equation 0,0]. In the present article, 
it will be shown that this can be done alternatively (but 
in a somewhat less strict way) using the gauge-invariant 
form of the transverse approach. For this purpose, we 
consider a quasi-one-dimensional (QID) two-band model 
with the impurity scattering taken into consideration. 

The two bands are taken to result from the site-energy 
dimerization in the highly conducting direction, and, con- 
sequently, the Bloch functions and all relevant vertex 
functions can be determined analytically (Sec. III). Using 
the equation of motion approach (which is found to be 
the generalization of the force-force correlation function 
approach), the intra- and interband optical conductivity 
are calculated (Sec. IV). In the intraband channel of the 
transverse correlation functions, the most singular pro- 
cesses in powers of l/u> n are collected, resulting in the 
optical conductivity which matches up the DC conduc- 
tivity obtained by the Boltzmann equations |lfj| or by 
the Landau response theory The resulting inter- 

band conductivity in the CDW ordered state is found to 
be consistent with the experimental observation. The- 
orywise, the semiconducting current-current correlation 
function 0, 0, Q is replaced by a slightly modified func- 
tion containing an additional factor which comes from 
the gauge invariant treatment of the diamagnetic cur- 
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rent contributions. Finally (Sec. V), the optical and DC 
conductivity are determined for a few typical Q1D cases. 



II. TRANSVERSE MULTIBAND RESPONSE 
THEORY 

The optical conductivity tensor a aa {uj) is a measure 
of the absorption rate for the (transverse) electromag- 
netic waves traveling across the crystal, and the measured 
spectra, together with the DC conductivity data and 
other transport coefficients, are an extremely valuable 
source of the information about the electronic subsys- 
tem. Although some aspects of the microscopic response 
theory can be found in the textbooks 0, IS 01 013 El El > 
there is no systematic microscopic solution to the multi- 
band optical conductivity problem. Actually, it is easy 
to determine macroscopic symmetry features of a aa {uj), 
even in a general case. In this respect, we shall combine 
the macroscopic symmetry features with the microscopic 
description of the electron-photon coupling functions (de- 
termined for a simple, exactly solvable Q1D electronic 
model) to develop a consistent microscopic multiband re- 
sponse theory. 



A. Optical conductivity tensor 

The optical conductivity analysis starts with the 
Hamiltonian El 

H = H§ M + H$ + HF t +HF t , (1) 

which comprises the bare photon term i?Q Cld , the bare 
electronic Hamiltonian Hq 1 and the first-order and the 
second-order electron-photon coupling term, Hf xt and 
i?2 Xt - The bare photon contribution is 

#0° ld = \ E t P q A + <0QUQ^} ■ (2) 

Here q and a are the wave vector and the polarization of 
the photon field Q qQ , P qct is the field conjugate to Q qQ , 
and Wqo = cq is the bare photon dispersion. The struc- 
ture of a typical Q1D tight-binding electronic Hamilto- 
nian H cl = Hq + Hi xt + is determined below. How- 
ever, notice that general symmetry properties of a aa (u)) 
discussed here do not depend on details of H el . 

To obtain a aa {uj), the retarded photon Green function 
is required Il4l | 

((Qqa'j Q-qa))t = -i0(t) < [Q qa (*), Q-q(O)]) (3) 

= e"" 4 — J d W e 1 " t ((Q qQ ;Q_ qQ )) w , 

with Q q „(t) = e^Q^e-^ and Q qQ (0) = Q qct 
representing the photon fields in the Heisenberg picture 



q,a q,a q,a Q q,a q,a q,a 

FIG. 1: The Dyson equation for the photon Green function. 
The bare second-order (first-order) electron-photon coupling 
leads to the second (third) term on the right-hand side, giving 
rise to the diamagnetic (current-current) contribution to the 
photon self-energy. 

at time t and t = 0, respectively. Using the equation of 
motion formalism, we get 

[u(u) + irj) -u; qce ]((Q qa ,;<2_ qQ )) w = H, (4) 

with the adiabatic term 77 — > 0" 1 " . The renormalized pho- 
ton frequency w qa is given by El 

w q Q = <4o+ fi dia,a + 47rI W w )> ( 5 ) 

where f^ia a ana - 47rn QQ (w) are, respectively, the dia- 
magnetic and current-current contributions to the pho- 
ton self-energy, shown in Fig. 1 (for the Q1D model under 
consideration, the explicit form of H aa (uj) and f^ia a ^ s 
given in Sees. IVB,C and Ref. EIj respectively). Com- 
bining the Maxwell equations with Eq. (J3J, it can be 
shown that the transverse dielectric function e a (ui) sat- 
isfies the relation 

[uj(uj + irj)e a (oj) - w q0 ] ((Q qQ ; Q-q a ))u> = K (6) 

with 

47ri 

e a (u)) = H a aa {uj). (7) 

UJ 

The optical conductivity defined by Eq. Q is 

«7oa(w) = (\.\ [Qdi»,a + 47rIIaaM], ( 8 ) 
4n[Ld + 17/) 

irrespective of whether the electromagnetic field is 
treated as a classical or a quantum field. This is a 
quite general and in many respects very useful result. 
In the general case, with several valence bands or with 
several scattering channels in a single band, f^ia a m " 
eludes one or more diamagnetic contributions (depend- 
ing on the number of bands intersecting the Fermi level), 
and TL aa (u)) represents all intra- and interband current- 
current correlation functions. 

The following general properties of the expression (JSJ 
are important for interpreting the measured spectra. 
First, there are at least two distinct structures in the 
optical conductivity spectrum Ke{a aa (uj)}, the first one 
is a delta function at u> = 0, related to the diamagnetic 
current, and the second one represents various contri- 
butions, including the exciton contributions if the short- 
range dipole-dipole interactions are present |3. IT5IT7Lll^ . 
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This is easily seen from 

Re{a QQ M} = i[n3 iBia +47rRe{II aa (0)}]*(w) 

— Im{n aa (w)}. (9) 

Obviously, in the normal metallic or insulating state the 
delta-function term vanishes [l3|. Therefore, any consis- 
tent treatment of the electron-photon coupling functions 
has to fulfill the relation 

nL, Q + 47rite{n QQ (o)} = 0. (io) 

The total optical conductivity in the normal metallic or 
insulating state can be written then in the form 

^aM = -[n Qa (w)-Re{n aQ (o)}]. (11) 

V 

Noteworthy, in the superconducting state the sum 

^["ka+47rRe{n QQ (0)}] (12) 

measures the weight of the missing area in optical con- 
ductivity spectra 0, |2(], while the single-particle 
contributions are still given by Eq. (|11|) . Finally, no- 
tice that in absence of local dipolar excitations (the case 
of the site-energy dimerization discussed in this article) 
the total current-current correlation function is the sum 
of only two contributions describing, respectively, the 
creation of the "free" intra- and interband electron-hole 
pairs, while the processes associated with excitons (i.e. 
the quasiparticles representing the bound electron-hole 
pairs) are not present. 

The optical conductivity determined by Eq. I|ll|l , with 
the current-current correlation function calculated by 
the equation of motion approach, is in the focus of the 
present analysis. The causality requirement 0, El (i- e - 
the Kramers-Kronig relations), the effective mass the- 
orem 13. p2l [23| and the gauge-invariance requirement 
0, ITU Il5| will be used to test the obtained results. A 
particular care will be devoted to the non-physical singu- 
larity at uj — related to the prefactor of oj^ 1 in Eq. (|llfl . 
and to the construction of the optical conductivity model 
with the correct behaviour in the to — > limit, giving rise 
to a unified description of the optical and transport phe- 
nomena. 



III. ELECTRONIC HAMILTONIAN WITH TWO 
QUALITATIVELY DIFFERENT SCATTERING 
CHANNELS 

Although our model is Q1D, the present response the- 
ory is quite general and the electronic Hamiltonian could 
represent an arbitrary multiband model. We assume 
that in addition to the bare electronic Hamiltonian (de- 
noted below by H ) there is a static single-particle po- 
tential (Hq) characterized by a commensurate wave vec- 
tor, which describes, for example, the scattering of elec- 
trons on the site-energy dimerization potential. The 



other single-electron scattering processes (on impurities, 
phonons, etc.) are represented by H[. For most of 
the questions discussed here the main effects of the two- 
electron interactions are taken satisfactorily into account 
through the effective mean fields included in H' or H[, 
|24j resulting finally in Hfi 1 = Hq + H' a + H[ in Eq. JTJ. 
jjext _ ffoxt^j^cxt C0U pi es the valence electrons to trans- 
verse electromagnetic fields. The resulting total elec- 
tronic Hamiltonian is 

H cl = H^ l + H cxt . (13) 

We start by diagonalizing the Hamiltonian Ho+Hq. As 
mentioned above, we consider the simplest Q1D model, 
where H' represents the site-energy dimerization in the 
highly conducting direction and in H[ only the impurity 
scattering is taken into account. Next, we determine the 
related electron-photon coupling functions. At the end 
of this section, the multiband current-current correlation 
function is introduced. 



A. Bare Hamiltonian 

The single-particle properties of the Q1D site-energy- 
dimerization model come from the exact diagonalization 
of the Hamiltonian 0] 

Ho + H' Q = [Mk)4 a c kCT + £c(k)cj cCT c kCT 

+A(c[ CT c kCT + 4 (T c kCT )]. (14) 

The bare electron dispersions of two subbands, artificially 
dimcrized along the highly conducting direction a, are 

£c. c (k) = ±2t a cos k • a — 2tb cos k • b, (15) 

with e c (k) = e c (k ± 7r/ax), and with the wave vec- 
tor k restricted to the new (reduced) Brillouin zone, 
— 0.57r/a < k x < 0.57r/a, — ir/b < k y < ir/b. t a and 
tb {ta 3> tb > 0) are the bond energies in the direction a 
and the perpendicular direction b, respectively, and A is 
the magnitude of the dimerization potential in the direc- 
tion a. Note that such potential corresponds to imperfect 
nesting in Eq. I|15fl in contrast to dimerization in all di- 
rections; the former nesting is chosen because it is more 
interesting in context of the conductivity studies. 

The transformations between the unperturbed states, 
the band index I — c,c, and the Bloch states, the band 
index L — C,C_, of the form 

L = Ecm^L, (i6) 

L 

lead to 

#o = ^^(k)4 CT L kff , (17) 

Lka 
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FIG. 2: The electron dispersions Ec,c(k x , 0.5n/b) in the re- 
duced (a) and extended (b) zone representation, for 2t a = 0.25 
eV and A = 0.02 eV. In Fig. (b), the solid and dashed lines 
correspond, respectively, to u(k) and v(k) in the relations 

CL = «(k)4. - "(k)cL £L = «(k)4 CT + «(k)^. The 
dot-dashed line is the Fermi level n in a typical metallic case. 

with the dispersions 



%,c(k) = i[ % (k)+e c (k)]±^ c (k) + A2, 
e £C (k) = e £ (k)-e c (k). (18) 

The transformation-matrix elements are given by 

(7k(c,C) f/ k (c,C) \ _ / u(k) «(k) 



U k (c,C) U k {c,C) 



where u(k) = cos ^ ^ , u(k) 



-w(k) u(k) 



, (19) 



2 . 2 
iary phase y(k) defined in the usual way 



sin — - — , with the auxil- 



tanyj(k) 



2A 

£cc(k) ' 



(20) 



The bands are shown in Fig. 2. Hereafter, the lower 
(conduction) band is assumed to be partially filled and 
the upper (valence) band is empty. For the lower band 
completely filled and A ^> tj,, the band structure corre- 
sponds to the commensurate CDW system, otherwise we 
have the metallic behaviour. 



B. Electron-photon coupling Hamiltonian 

Using the generalized minimal substitution method for 
the tight-binding electrons 0, we obtain that the 
conduction electrons described by the Hamiltonian l|14|) 
are coupled to the external electromagnetic fields through 

HCXt = EEE ffl o(MC^- ( 21 ) 

/ ker qct 

where, to the second order in the vector potential A a (q), 
o T w„ x dedk) e , , . 19 2 e;(k)/ e\ 2 ,,, . 

(22) 

The photon annihilation operator A qa enters in Eq. i|22|) 
through 



A*(q) 



4ttc 2 

5> Q (q-q')A Q (q'), 
q' 



where 



(23) 



is the electromagnetic displacement field of Eq. J2J) _11]. 
Finally, in the Bloch representation the coupling Hamil- 
tonian becomes 



H 



r' + fff = --V4»(q)i (-q) 

r — ' 



2mc 2 



X)^(q)7aa(-q;2), 



with 



^a(q) 
7aa(q; 2 ) 



LL' ker 



(24) 

(25) 
(25') 



LL' ko 



representing, respectively, the current density and Ra- 
man |l5l |22| density operators. The structure of the re- 
lated vertex functions J^ L (k) and 7^ (k; 2) is given in 
Appendix A. 



C. Single-particle scattering processes 

We consider here only the intraband impurity scatter- 
ing processes in the perturbation 

L kk'o- 

for which one usually assumes V LL (k, k') = V LL (\s.— k') 
0, . This form of H[ is consistent with the regime in 
which e cc (kp) 3> H/t (t is the relaxation time defined 
below). The generalization is straightforward. 
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FIG. 4: The generalized correlation function in notations of 
Fig. 3. 
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FIG. 3: (a) The Bloch representation of the current-current 
correlation function of the two-band model 11141 1 . with the 
current vertices given by Eqs. 1A21 and (A2'). The shaded 
square is the electron-hole propagator defined by Eq. 1281 . (b) 
The intraband current-current correlation function at high- 
frequencies (the leading term in powers of (H[) 2 /to), (c) The 
structure of the high-frequencies self-energy term. The direct 
processes are associated with the creation of an electron-hole 
pair with wave vectors k and k + q (q is the wave vector 
of the external field), while the indirect processes involve the 
impurity-assisted electron-hole pairs with wave vectors k and 
k' + q ~ k', or k + q ~ k and k' (k' — k is the momentum 
relaxed on impurities, and the crosses represent the impurity 
scattering) . 



D. Current-current correlation function 

In the equation of motion approach [Tlj used in the 
next section, the starting point is the current-current cor- 
relation function of Eqs. JSJl-dJ shown in Fig. 3. It is 
defined by 0,110 d 



n« Q (q,<) = — ((J Q (q); J Q (-q))) t 

= -^ 7 e(i)([j Q (q,t),J a (-q,0)]). (27) 



Here the current operator J a (q) includes all intra- and 
interband current density fluctuations, as seen from 
Eq. (J22J- According to Fig. 3, II QQ (q, t) comprises two 
intraband and two interband contributions, and the prob- 
lem is reduced, as will be seen immediately below, to 
the self-consistent calculation of intra- and interband 
electron-hole propagator in presence of the perturbation 
H[. In Eq. as well as in Eqs. QgJ, (JSOJ), (33) 

and H34fl . the usual notation for the retarded correlation 
functions is used: ((A,B)) t = -iQ(t){ [A(t), B(0)] ), with 
A(t) being the operator A in the Heisenberg picture. 



IV. EQUATION OF MOTION APPROACH 
A. Generalized correlation functions 

The retarded electron-hole propagator 

pf L '(k,k + ,kV,k',<) = <<4X +q ^t +qCT £kv)) t 

ee -ie(t)([(4 CT L' k+qCT ) t , (Lt, +qff Lw) ]) 
1 f°° 

= e-"*— y dwe^ £ (k,k' + ,k+,k»(28) 

(k + is the abbreviation for k + q) is the central quantity 
to all long-wavelength correlation functions, as can be 
seen from Fig. 4, or from the expression 



Xf„ 



(q> f ) = w S £ fLL '^ k + q)5 L,i ( k ' + q. k ') 



LL' kk'o- 
■\LL' 



xV[ L (k,k +) k' + ,k',i), 



(29) 



which represents a generalized long-wavelength correla- 
tion function, with f LL '(k, k + q) and g L ' L {k' + q,k') 
being the charge, current, Raman or some other ver- 
tex functions. In Eq. H28fl the abbreviation (AB) t = 
A(t)B(t) is used. 

The symmetry of the f LL '(k, k+q) and g L ' L (k'+q, k') 
vertices, together with the nature of the singularity 
of the leading term in the perturbation series, deter- 
mines the summation rule for the related Feynman di- 
agrams. The simplest way to collect the most singular 
diagrams is to consider the equations of motion connect- 
ing T>\ L (k, k + , k' + , k', t) with the correlation function 



T>% L '(k,k + ,k',,k',t) defined 



V% L '(k,k + ,k' + ,k',t) 



(30) 



The direct calculation gives the exact relation 

h[V^ L \k, k + ,u)] - X V\ L ' (k, k+, k' + , k', w) (31) 
= M k , k , [/ £ (k) - /z/(k+)] + P 2 Li '(k, k + ,k' + ,k',u;). 

Here 



h\v^ L '(k,k + ^) 



h(u + irj) + E L (k) - E L >(k + ) 



(32) 



is a useful abbreviation. T>% L (k, k + , k' + , k', to) is the 
Fourier transform of T>\ L '(k, k + , k+, k', t), and /i(k) = 
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/(^(k)) = {I^L^l = [1 + exp{(3[E L (k) - M ]}] 1 is 
the Fermi-Dirac function, with (3 = l/{knT). 

The way to evaluate L (k, k+, k' + , k', t) depends on 
the choice of representation of this function. There are 
two alternative ways, 

P 2 Li '(k,k+,k' + ,k',i) (33) 

= -ie(t)<[[4x +qCT)J ff'] t , {pl +vt L kfa )^ ) 

or 

P 2 Li '(k,k + ,k;,k',i) (34) 

= -ie (t) ( [ [4X+ q - o . (4 +q ,ikv) _ t ] ) , 

leading to two different self-consistent schemes (one de- 
scribed below and another encountered in the longitudi- 
nal response theory 0), both giving, as will be argued 
below, the same result. Again, [A, B] t is the abbreviation 
for [A(t),B(t)}. 

It is important to realize at the outset that the first 
term on the right-hand side of Eq. (|31|l is significant for 
all interband correlation functions x/, 9 (q, in Eq. 
irrespective of the vertex symmetries. The probability 
for the direct creation of the interband electron-hole pair 
is proportional here to / L (k) -/ L /(k+) w f L (k) - f L > (k) , 
so that all occupied states in the conduction band(s) 
are equally important. For vertices j LL (k, k + q) and 
g L L {k' + q, k') taken to represent the intraband current 
vertices, the correlation function x/,g(q>*) in Eq. 
becomes the intraband current-current correlation func- 
tion n™* ra (q, t). According to the discussion of Sec. II A, 
the related intraband optical conductivity is ruled by 
the prefactor lo^ 1 in Eq. (|lip. The direct processes in 
Eq. H31|) (see the first term on the right-hand side of 
Fig. 3(c)) related to /c(k) — /c(k+) — > are insignifi- 
cant, and IP c j^ ra (q, t) can be adequately described by the 
second, indirect term in this equation (the second term 
in Fig. 3(c); see also Fig. 5), with X>f- L '(k,k + ,k' + ,k',i) 
given by Eq. . 

Before turning to the evaluation of II^* ra (q, t) it is in- 
teresting to contrast this conclusion with its analog in 
the longitudinal response theory. In the longitudinal ap- 
proach, the long-range Coulomb forces are activated, and 
the first term in Eq. H31J1 dominates the intraband charge- 
charge correlation function, even in the dynamic limit, 
since the small probability for the intraband electron-hole 
pair creation, proportional to /c(k) — /c(k+), cancels 
out the q~ 2 singularity of the long-range forces, which 
is the well-known RPA result. The indirect scattering 
processes, on the other hand, are proportional to the ef- 
fective intraband charge vertex, analogous to the effec- 
tive intraband current vertex introduced below, Eq. 1)41 [I. 
Since the bare long-wavelength intraband charge ver- 
tex satisfies q(k + q, k) ea q(k, k) = e, where e is the 
bare electron charge, the effective intraband charge ver- 
tex for the indirect scattering processes vanishes, because 
q(k, k) — q(k' , k') = 0. Since the longitudinal and trans- 
verse approaches are to be equivalent, this means that 
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FIG. 5: The direct intraband (1 — ► 2), indirect intraband 
(1 — > 3 — ► 4), direct interband (1 — > 5) and indirect interband 
(1 — > 6 — > 7) optical excitations. In this article the second and 
third processes are considered, the indirect intraband, using 
the exact summation of the impurity scattering processes, and 
the direct interband, using the phenomenological treatment 
of the impurity scattering. The dotted line is the photon 
dispersion, and the dashed line corresponds to the impurity 
scattering. 

the direct longitudinal processes are to be equivalent to 
the indirect transverse processes in the intraband chan- 
nel, while the contributions of both the indirect longitu- 
dinal and the direct transverse scattering processes are 
to be negligible. This issue will be further discussed in 
Sec. IV B 5. 

The above transverse approach can now be applied to 
the current-current correlation function of the two band 
model, rewritten in the form 

lWq,t) = ^EE ^ L '(k)jf £ (k') 

LL' kk'o- 

xPf L '(k,k+,k;,k',t). (35) 

Here the vertices f iL '(k, k + q) and g L L {k' + q, k') in 
Eq. JUJ) are replaced by J^ L '(k,k + q) « J £i '(k) and 
J^' L (k' + q,k') « J^' L (k'), respectively. Optical pro- 
cesses relevant to the two-band model (including the in- 
direct interband ones not considered in the present anal- 
ysis) are illustrated in Fig. 5. 

In order to make presentation of the results more trans- 
parent, we shall first determine the intraband contribu- 
tions, and then give the analysis of the interband optical 
excitations. For the sake of brevity, in the next section 
the (intra)band index C will be omitted (£?c(k) — > E(k), 
j£ c (k) J Q (k), with Cl a -> 4J- It should 
be noticed that the results obtained below for the in- 
traband optical conductivity are quite general, i.e. they 
cover various physically different regimes. As explained 
in Sec. IV B 2, when the electron filling of the dimerized 
band varies between and 1, the electronic system trans- 
forms from an electron-like semiconducting, through a 
metallic, into a hole-like semiconducting regime. A more 
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detailed discussion of this issue is given in Sec. V, in con- 
text of the total optical conductivity. 

B. Intraband optical conductivity 

According to the aforementioned arguments, the intra- 
band optical processes are described by the equations 

hV^ 1 (k, k+ , w)2?i (k, k+ , k' + , k' , u) 

= V 2 (k,k + ,k' + ,k',uj), (36) 
?i© - 1 (k',kV ) a;)© 2 (k,k + ,k' +) k / ,a;) 
= -P 3 (k,k+,k' + ,k» 

+ h i [ [ c L c k+ q <r , H'] , 4, +qo .c k v] } , (37) 

where T>Q 1 (k', k' + , uS) is the intraband term in Eq. I|32|) 
and 2?3(k, k+, k^, k',o;) is the Fourier transform of the 
force-force correlation function 0, 0, |2f| 

£> 3 (k,k + ,kV,k',t) (38) 
= -ie(t)([[cl a c k+ ^,H'] Q , [4, +qff c k v, H']_ t }). 

The second term in Eq. (|3T)l is the ground-state average 
of the four-operator product at t = 0. It is off-diagonal in 
the Bloch representation, and its value can be obtained 
by putting E(k' + ) - E(k') rj in the left-hand side of 
Eq. I|37|) and then taking the formal limit to — > 0. The 
result is 

H [ [ct a c k +q*, H'] , 4, +qCT c kv ] ) 

w2? 3 (k,k + ,k^,k / ,w = 0). (39) 

For the impurity scattering, we have the cancellation of 
this constant term with its counterpart obtained by the 
electron f± hole replacement (see Eq. (14711 . for exam- 
ple), but for a time-dependent perturbation H[, Eq. ljoT?|) 
achieves an interesting structure, as pointed out in Refs. 

Due to the large velocity of light, the energy and wave 
vector transfers in the external points of the intraband 
current-current correlation function fulfill fku 3> E(k',) — 
E(k'), and the factors T>o(k, k + , uj) and T>o(k' , k' + , uS) in 
Eqs. <|36|) and (|37ll . representing the propagator of the 
virtual electron-hole pairs (related to the process 1 — > 3 in 
Fig. 5), can be replaced by 1/oj. The resulting intraband 
correlation function becomes 

^ ' kk"kiki'o- 

x [j a (k'{) - J a (ki)] (v(k - k")V(k'( - kx) 

x [X> x (k, k", k'/, ki, w) - 2?i(k, k", k'{, ki, 0)] ^. (40) 

(. . .) is here and subsequently the average over the im- 
purity sites 0. 




FIG. 6: (a) The indirect contributions to the intraband 
current-current correlation function in the equation of mo- 
tion approach. The open square is the effective current ver- 
tex, and the dashed line represents the related force-force 
correlation functions (see, for example, the discussion of 
the force-force correlation function method in Sec. 8.1.B of 
Ref. 13 )• (b) The leading term in the high-frequency limit. 
On the right-hand side, the products V(k — k")J a (k)/(hw) 
and V(k — k")J a (k")/(hu) are encircled into the effective 
vertex. The dot is the bare intraband current vertex. 

The diagrammatic representation of II" 1 * ra (q, u>) is 
given in Fig. 6(a), with the structure of the leading term, 
proportional to (iJ() 2 /w, shown explicitly in Fig. 6(b). 
The expression 

j a (k,k") = V(k ? ~ k " ) [-JafrWaQO] (41) 

can be recognized as an effective current vertex for the in- 
direct intraband photon absorption/emission processes, 
and is a sum of two terms encircled in Fig. 6(b). Notice 
also that Di(k, k", k'{, k x , u) - X>i(k, k", k'{, ki, 0) is pro- 
portional to lj, canceling out the factor of in one of 
the effective vertices. 

It is important to notice that in the leading term 
(Fig. 6(b)) the expression (j4"U|) is identical to the result 
of the force-force correlation function method • While 
the latter method is usually limited to the examination 
of this leading term, or to the summation of irrelevant 
higher order diagrams which results in the well-known 
Hopficld formula |JJ, the present approach is focused in- 
stead on the exact summation of the most singular con- 
tributions in powers of (H[) 2 /uj and should be regarded 
as a generalization of the standard force-force correlation 
function approach. 



1. Proper electron-hole representation 

After determining the structure of the external points 
(i.e. the effective vertices) in the diagram shown in 
Fig. 6(a), we have to find the internal structure of the dia- 
gram. The latter is represented by the indirect (impurity- 
assisted) electron- hole propagator £>i(k, k", k", ki, w), 
characterized by the momentum transfer k — k", rather 
than by the negligibly small external momentum transfer 
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(a) 



two equations connecting the electron-hole propagator 
Pi(k,k",ki',ki,w) with the kernel F a {k'{, ki, w). The 
solution is based on the expansion 




7r 

(d) k, 



FIG. 7: (a) A few leading contributions included in the self- 
consistent treatment of Eq. (1421 . (b) Typical contributions 
omitted in the self-consistent scheme, (c) The self-consistent 
equation of the kernel ^(k", ki,w) and (d) of the electron- 
hole propagator T>\ (k, k" , k" , ki , uS) . The diamond represents 
the electron-hole self-energy consisting of the single-particle 
self-energy and vertex-correction contributions. 



I) 1 (k 1 k" 1 k' 1 ',k 1 , W )=^P! 2n) (k,k",k' 1 ',k 1 ,,), (45) 

n=0 

where T)^ n \k, k", k", ki,o;) includes only the self- 
consistent terms proportional to (H[) 2n / ui n . 

As easily seen in the longitudinal analysis, the self- 
consistent expression for 2?i(k, k", k", ki, u), with k" — > 
k + q, plays the leading role. Accordingly to com- 
pare both response theories, one needs the self-consistent 
scheme for ^.(k", ki, ui) and the recurrence relations for 

Pf n) (k,k",ki / ,k 1 ,cc;) (see Sec. IV B 3). 

For pedagogical reasons, it is convenient first to con- 
sider the zeroth order contribution to l|43|) and define the 
effective number of conduction electrons and the related 
electron-hole damping energy. 



2. Htgh-frequency limit 



k + — k = q of the direct processes in Eq. (|31() . Its in- 
ternal structure is determined here by the self-consistent 
solution of the exact equation 

hV^\k, k", w)£>i(k, k", k'{, ki, w) (42) 
= M k , k A»,k»[/( k ) ^ Z( k ")] +P 2 (k,k",k' 1 ',k 1 ,^). 

2?i(k, k", k", ki, uo) represents the electron- hole pair cre- 
ated by the indirect photon absorption of Fig. 5, which 
means that the energy transfer hio is close to the electron- 
hole pair energy E(k) — E(k"). The wave vectors k and 
k" (ki and k", as well) are independent of each other, 
and therefore the first term in Eq. (|42|l dominates the 
low-energy physics, preferring the optical processes be- 
tween the states E(k) w fi and E{k") w \i. 

In combining Eq. I|42(l with the equation of motion for 
2?2(k, k", k", ki, u) (see Eq. l|53|) for more detail) it is suf- 
ficient to collect only the terms in the summation which 
are relevant to the self-consistent equation for the kernel 
J- a (k!{, ki,w) in the current-current correlation function 
(see Figs. 7(c,d), and the criterion of the validity shown 
in Figs. 7(a,b)) 

n ^ ra(w) = ~whf E ([JM)-J«(M] 



: [^.(ki'.ki.w) -^.(kjf.ki.O)] Y (43) 



where 



^(k'/,^^) = £V(k - k")V(k'( - k x ) 



kk" 



x[j a (k) - J a (k")]Pi(k,k",k; , ,k 1 ,w). (44) 

The idea of the present approach is the self-consistent 
treatment of Eqs. (|42|l and (|44|) . which represent 



The direct calculation gives the first term in the ex- 
pansion l|45|) 



vf\k,k",KM^) 



(46) 



= <5 k , kl <W/[.f( k ) - /(k")]Z>o(k,k'». 

The related contribution to the current-current cor- 
relation function (which is a good approximation for 
n'aa ra (w) at high frequencies) is given by 

1 1 



n 



intra, (0) 



( Ja(k)(j„(k)-J a (k")) 



kk" 



x[/(k)-/(k")](ink-k")| 2 ) 

xi[D (k,k /, ,w)-X>o(k",k ) «)] (47) 

w 7^17 E ^(k)(J«(k) - J a (k")) 
(Tuj) v kk „ CT 

x [f(E(k)) f(E(k) + M] (\V(k - k")| 2 ) 
1 



x - [X> (k, k", u) + D (k", k, w)] . (47') 

For the impurity scattering processes, the real part of this 
function is negligible, while the imaginary part is given 
by 

Im{n^°>H} « -77^72 £ E J '( k ) ( 48 ) 



kcr 



x [f(E(k)) - /(£(k) + fia;)] 
where the electron-hole damping energy is 



r(k,w)- 



^ = X;(|y(k-k")| 2 )(i-^g 



r(k,cj) 



x(-)-Im{2> (k,k'»}. (49) 
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a 

> 



0.8 



0.6 



0.4 



0.2 




electron- like 



0.5 1 

5 (electron doping) 

FIG. 8: The effective number of conduction electrons 15111 
shown in the dimensionless form, for 2t a = 0.25 eV, tb — > 0, 
T = 10 K and A = (dotted line), 10 meV (solid line). The 
dot-dashed line is the dimensionless concentration of conduc- 
tion electrons Von = 2ak-e/n = S, i.e. the prediction of the 
free-electron approximation. Vo is the primitive cell volume 
and m aa = h 2 /(2t a a 2 ) is the mass scale parameter. The dop- 
ing regions with the electron-like (9nf^ ra a /d5 > 0) and hole- 
like (dnfntra, a /d5 < 0) conductivity are also indicated. Inset: 
The effective number (m aa /m)Vonf* tra a for 2t a = 0.25 eV, 
t b -> 0, T = 10 K and A = 0.25 eV. The dot-dashed line is 
the prediction of the free-hole approximation, 1 — Von. 



Furthermore, in this case, the frequency dependent part 
in l/r(k, u>) is negligibly small, and we can put the 
average over the Fermi surface (l/r(k, 0)}fs = 1/ T m 
Eq. ggj instead of l/r(k,w). We finally get 



2„cS 



im{ir Q n r a < (0) M} 



en 



intra. a 



with 



n 



off 



dm 

dE(k) 



(50) 



(51) 



being the effective number of conduction electrons. The 
behaviour of nf^ tla a with band filling S is shown in Fig. 8 
for a few typical values of the ratio A/(2i a ), and com- 
pared to the prediction of the free-electron and free-hole 
approximation. Notice that for A > 2t a » 2tb one ob- 
tains the well-known result n^ tm a cx sin k f x 2a — sin Sir 
(2a is the primitive-cell parameter of the dimerized lat- 
tice). 

Beyond this approximation the expression (|47|l re- 
quires the evaluation of two coupled integrations over k 
and k". In the low-dimensional electronic systems, where 
the van Hove singularities in the band structure may play 
important role, this is not a trivial task. 

The most outstanding advantage of the present ap- 
proach is the fact that the effective current vertex 
j Q (k, k") consists of two terms, which, together with 
two terms in j a {k'l,k\), give four contributions to 



j-jmtra,(o)^^ twQ Q £ tnem representing the so-called self- 
energy contributions, and the other two the vertex cor- 
rections (see Fig. 6(b)) 0,E3- We shall show next that in 
the nth order term in (|45|) the single-particle vertex cor- 
rections and the single-particle self-energy contributions 
are also treated on equal footing, even if the relaxation 
processes on impurities are treated beyond the approxi- 
mation l/r(k, cj) » 1/t. 



3. Kernel in the low-frequency limit 

The effects of the impurity scattering on the correlation 
function l|37|) can be exhibited in the following way: 

V 2 (k, k", k'{ , k x , w) = £ V (q) [£>r (k, k" - q, k'{ , kr , u) 



-X> 1 (k + q,k",ki , ,k 1 ,a;)]. 



(52) 



When combined with this expression, Eqs. (|42|1 and H4bj) 
lead to 



Di(k, k", k'/, kx, u) « V { ° ] (k, k", k'{, ki, u 



(53) 



ipo(k,k", W )^|^(q)| 2 



X - {V (k, k" + q, u) + V (k + q, k" , u) } 

x [X>i(k, k", k'{, k l5 u) - X>i(k + q, k" + q, k'/, ki, u)] . 

On the right-hand side of this equation, only the self- 
consistent terms are taken into account. The first term 
in the brackets represents the single-particle self-energy 
contributions and the second one the single-particle ver- 
tex corrections. It is important to remember again that 
these corrections are the largest for E(k) w E(k"), so 
that £>o(k, k", u) can be approximated by w" 1 . As a con- 
sequence, the solution of this equation can be sought in 
powers of \V (q)\ 2 / cu . However, due to the dependence on 
k + q and k" + q of the vertex-corrections term, we have 
to turn back to the kernel (|44|) and apply several changes 
to its dummy variables to obtain the self-consistent de- 
scription of T a (k![, ki, to) and the desired recurrence re- 
lations between the contributions V^ n ^ (k, k", k", ki, u>). 
The kernel is described by 

7(0)1 



= ±Y\v(k-k>iy 2 



]T|^(k-k")| z [J Q (k")- j a (k)] 



kk" 



xP 1 (k,k",k' 1 ',k 1 , W )S(k,k", W 

where 

t£E(k, k",u) 

= -E(i^- k ")i 2 )f 1 - 



(54) 



(55) 



D Q (k,q,o;) 



E 



(\V(q-k)\ 2 ) 1- 



J«(k" 

Jjk))h Vo ^ k ^ 
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is the electron- hole self-energy and T^(k!{, ki, uj) quirement 
is the kernel corresponding to the replacement of 

(0), 



2?i(k,k",ki / ,ki,w) in Eq. O by V?> (k, k", k'{, ki, uj). 
The only approximation made in the derivation of Eq. 
(55) is 



J a (k - k" + q) - J Q (q) « [j Q (k) - J Q (k")] 



Jq(q) 

Ja(k")' 

(56) 



which allows a simple description of the vertex correc- 
tions in ,F a (k", k 1; uj), and which treats correctly the 
disappearance of the forward scattering contributions 
(k w k") to both £(k,k",w) and T a Qx!l, k x , w). As 
briefly discussed at the end of Sec. IV B 5, this approx- 
imation is directly related to the restrictions enforced by 
the the continuity equation for the charge density. Simi- 
larly, the definition l|45|) , together with the self-consistent 
relations (|53|l and (|54|l . gives rise to the recurrence rela- 
tions illustrated in Fig. 7(d): 



(57) 



-E(k,k'» g(2ra _ 2 



-£(k,k'» 



'(kk",^',^) 



Re{ar a (-^)} 
Im{<-(-^)} 



Re{4 ntri ' 
-Im{< ra H}, 



,_intra / 



(60) 



provided that r(w) = t. If t(uj) is frequency dependent, 
the corresponding Re{E(k, k", uj)} is non-zero, but the 
result (|59J) is still acceptable. For Re{E(k, k", u>)} not 
too large we can introduce the effects of Re{£(k, k", to)} 
through the mass redefinition m — > m{ui) through the 
Kramers-Kronig relations. The result is the generalized 
Drude formula for the intraband optical conductivity. 

Here we show two important results. First, the 
memory-function approximation, which in the traditional 
form has not been found to be transparent, can be un- 
derstood as a simple replacement E(k, uj) — > i/r(cj) of 
the exact-summation result given below. Second, the 
memory-function results are acceptable even in the cases 
with the pronounced optical excitations across a gap (or 
pseudogap) (where n'?^ tla Q <C n; see, for example, the 
hole-like semiconducting regime in Fig. 8 at S « 1), pro- 
vided that the electron group velocity v a (k) = J a (k) / e 
in Eq. 1)5 1H is determined using the relation 1A5(I . The in- 
traband conductivity spectrum obtained in this way is re- 
lated with the interband conductivity spectrum through 
the well-controlled conductivity sum rule |15| . 



4- Memory- function approximation 



Exact summation 



The simplest way to solve Eq. <|54|) is to replace 
E(k, k", uj) by its imaginary part averaged over the Fermi 
surface, i/t(uj). As mentioned above, for the impurity 
scattering processes, the real part of £(k, k",w) can be 
ignored. Even if Re{E(k, k", uj)} is not small, we can 
turn back to the electronic Hamiltonian and tray to in- 
clude the Re{E(k, k", u>)} effects into the effective single- 
particle Hamiltonian, and to diagonalize it, as we did here 
with the scattering processes on the dimerization poten- 
tial Hq. The real part of the new self-energy E(k, k", uj) 
is minimized in this way, with only the imaginary part 
playing an important role in Eq. 1)570. 

In this case, we obtain the expression 



UJ + \/T{ljj) 



, (58) 



which leads to the well-known results of the memory- 
function approximation 0, |2|| 



m uj + i/t(u>) 



- J, . y (59) 

uj m uj + i/t(uj) 

with H/t(u>) being the intraband memory (relaxation) 
function. This result is consistent with the causality re- 



in the case where the wave vector dependence of the 
imaginary part of E(k, k",u>) is significant and the real 
part of E(k, k",u>) is not too large the full recurrence 
relations for X>j 2r ^ (k, k", k'/, ki , uj) can be used to obtain 
the intraband current-current correlation function. The 
result is 



kk"o- 



7 E 400 



kk"<r 



hhuj + hz(k,k",uj) 
df(k) l 
dE{k) fao + hZ(k,k",uj) 



x(\V(k-k")\ 2 ) 1 



J a (k) 



x i [V (k, k", uj) + D (k", k, uj)] . (61') 



Physically the most important case corresponds to the 
approximation E(k, k",uj) « E(k, k, uj) — S(k, uj). In 
this case, we obtain 



nL n r a M 



ktr 



df(k) s(k,c) 

dE(k) w + E(k,w)' 



(61'' 
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The final form of the optical conductivity comes from 
Eqs. (61") and (11) 



i 1 



9/(k) 



(62) 



which is the result identical to the result of the longitu- 
dinal response theory. 

The longitudinal response theory gives a simpler way 
to obtain the same txj^ ra (a;). The difference between the 
two approaches is in the way how the continuity equa- 
tion connecting the charge density and current density 
fluctuations is treated. The consideration of the direct 
processes of the wave vector q in the longitudinal ap- 
proach allows a more precise treatment of the continuity 
equation, in the way analogous to the Landau response 
theory [ill Il2| . But here only an approximate solution 
is possible, since the theory is formulated in terms of the 
indirect intraband processes (of the wave vector k — k"). 



6. Zero- frequency limit 

First significant consequence of the present exact sum- 
mation method is the behaviour of the intraband optical 
conductivity in the zero- frequency limit. When the real 
part of E(k, u>) is small enough, we can write 



S(k, 



iS"(k,0)=i/r(k), 



(63) 



resulting in the DC conductivit y w hich is equal to the 
well-known Boltzmann result [?lllO| 

£"(0)^° = (-)I^j2( k )|ffl T ( k ) 



ktr 



e TO m -eff 

mVb m a 



intra, a ' 



(64) 



Here 



n cff 



aa \r „cff 



, ,m aa 2 3/(k) r(k) 



(51') 



is the effective number of conduction electrons shown in 
the dimensionless form and To = r(k = 0) is the tem- 
perature dependent k = relaxation time. The relation 
B62JI . together with Eq. (|66[) . gives the complete descrip- 
tion of the optical conductivity in a general multiband 
model, with the symmetry of the intra- and interband 
current vertices playing an essential role. Thus, Eq. i|64|) 
provides the direct link between the low-frequency op- 
tical conductivity and various transport coefficients not 
only in the single-band but also in the multiband models. 



C. Interband optical conductivity 

The approximation in which the 2?2(k, k", ki, oj) term 
in Eq. Q31[) is omitted leads to the ideal interband current- 
current correlation function and to the ideal interband 
conductivity characterized by a sharp threshold at the 
energy %(k F ) - E c (k F ) 0,0,11 EH- The former one is 
given by 



~ v h \^-%c(k)+i 
/c(k)-/ c (k) 



ftu + £fec(k) + i/w? 



ifvq 
(65) 



with i?cc(k) = i?c(k) - -Ec(k). The latter one comes 
from Eqs. ifTTj) and l|63|l 



winter/, ,\ 



i 1 >p frw\J§-° (k)\ 2 
wVf S cc (k) 

ktr — 



2[/ c (k)-/c(k)] 
hu + ihrj - E^ c (k)/(Huj) 



(66) 



with i] — put in the numerator and rj — > + in the de- 
nominator. The leading term in the interband electron- 
hole self-energy 



K^ll 1 (k, lu) = -J2 



(|^' L '(q-k)| 2 )il? "'(k,q,c) 



1 



+ <|V^(q-k)f)-2^ (q,k,c) 



(67) 



comes from the single-particle self-energy contributions. 
A reasonable generalization for the interband optical con- 
ductivity is given by Eq. I|66|l with the replacement r] by 
Im{E S c(k,w)}. 

The interplay between the self-energy and vertex- 
corrections terms in Ecc(k, lu), the correct treatment 
of the indirect interband optical excitations, as well as 
the role of the effective mass theorem in resolving all 
these issues will be explained elsewhere 01 • It should 
be noticed here that the Landau-like function (|66|l gives 
the in-gap optical conductivity slightly different from the 
corresponding Lindhard-like function, as easily seen by 
comparing Fig. 9 with Fig. 4 reported in Ref. fTs]. The 
understanding of the difference between these two results 
is of general importance as well, and will be discussed in 



detail in Ref. 



mpc 

m 



V. COMPARISON WITH EXPERIMENTS 

In the simples limit, £cc(k, oj) ~ iri tltra and 
Scc(k, lu) ~ iFintor, the optical conductivity of the 
present two-band model, 



a aa {u) 



intra 



(7 



inter / 
aa v 



LU 



(68) 
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A =0 

A = 2 meV; T = 64 K 
A = 4 meV; T = 56 K 
A = 6 meV; T = 42 K 
A = 8 meV; T = 24 K 
A=10meV;T = 5K 
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FIG. 9: The development of the normalized total single- 
particle optical conductivity with temperature in the ordered 
CDW state: A(T) = A y/l - T/T C dw, A = 10 meV, 
Tcdw = 66 K, fc F = 0.57r/a, 2t a = 0.25 eV, t b -> and 
fiFintra = fiFintcr = 5 meV. o\a labels the DC conductivity 
at A = 0. 



1 08 : io 



S 0.6 



0.4 



c 0.2 




T / T 

' CDW 

FIG. 10: The temperature dependence of the normalized 
effective number of conduction electrons in the ordered 
CDW state for different values of exponent j3: Ao = 10 
meV, Tcdw = 66 K, fc F = 0.57r/a, 2t a = 0.25 eV, 
tb —* and fiFintra = ftFinter = 5 meV. Notice that 
this figure also represents the temperature behaviour of the 
normalized DC conductivity, since a£f (T) / (?aa (Tcdw) ~ 

"Ltra.a ( T ) I "-hftra.a (TCDW ) • 



we contrast the interband conductivity found here to the 
oversimplified semiconducting optical conductivity usu- 
ally encountered in the textbooks pi 0, Q ■ 



A. DC conductivity in the ordered CDW state 

The temperature dependence in k^T, A(T), Fi n tra(T) 
and Tenter is responsible for the transfer with increas- 
ing/decreasing temperature of the optical conductivity 
spectra between the intraband and interband channels. 
This effect is particularly large in the vicinity of the 
metal-to-insulator phase transition. It should be noticed 
that, in the approximation in which r(k) — > tq = l/r; ntra 
is set in Eq. (51'), the temperature dependence of er°^ is 
given by the product of the effective number of conduc- 
tion electrons nf^ tIaa (T) and the relaxation time tq(T). 
We have two adjustable parameters at any tempera- 
ture, and the analysis of the DC conductivity data is 
thus possible only by the combination with the optical 
conductivity measurements. The latter method allows 
also the determination of the magnitude of CDW order 
parameter A and the critical exponent (3 in A(T) = 

A (l — T/T CDW ) , as well as the bond energy t a and the 
damping energy r int cr. 

In Figs. 9 and 10 the optical conductivity normalized 
to the DC conductivity at A = and the temperature 
dependence of the DC conductivity are shown for typical 
values of the parameters. 



B. Optical conductivity of a simple semiconductor 

The typical result for the interband conductivity in the 
ordered CDW state is shown in Fig. 11 (solid and dot- 
ted curves), and compared to the data measured in the 
blue bronze K0.3M0O3 (diamonds) [Bj]. Notice that the 
gauge- mvariancc factor huj/Egc(k) in Eq. JS3) ensures 
the disappearance of cr™* er (0) at T — » 0, independently 
of the value of the phenomenologically introduced damp- 
ing energy hTi nteT . The dashed curve is the prediction of 
the usual optical model for semiconductors 0, 0, 0] , with 
the factor Huj/EccO*) absent, which is characterized by 
a significant (but non-physical) contribution to at 



is a function of the Fermi wave vector kp, three band 
parameters, t a , tb and A, and two damping energies, 
^Tintra and fiFj n te r . For k-px ~ 0.57r/a, the model il- 
lustrates optical properties of various Q1D imperfectly 
nested CDW systems (including both the ordered CDW 
state and the pseudogap effects at temperatures above 
the critical temperature Tcdw)- In this section, we shall 
briefly discuss a few qualitative results important to the 
Q1D CDW systems. First, the temperature dependence 
of the DC and optical conductivity in the ordered CDW 
state is discussed for the strictly ID case (tb — > 0). Then 



VI. CONCLUSION 

In this article, the response of the conduction electrons 
in a Q1D multiband model to the transverse electromag- 
netic fields has been studied in presence of two scatter- 
ing mechanisms. In contrast to the coherent scattering 
on the site-energy-dimerization potential, the impurity 
scattering processes give negligibly small contribution to 
the real part of the electron self-energy, but dominate in 
the relaxation processes in both the intra- and interband 
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0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 
energy (eV) 

FIG. 11: The interband optical conductivity 16611 (hrj — > 
fir intcr ), for 2t a = 0.4 eV, t b -> 0, A = 0.09 eV, fe F = 0.57r/a, 
T = 10 K and Winter = 0.3 eV (solid curve) and 1 meV 
(dotted curve). The dashed curve is the prediction of the 
usual (gauge- non-invariant) optical conductivity for semicon- 
ductors, calculated using the same band parameters, with 
ftr intcr = 0.15 eV. <r° a c is the DC conductivity of the A = 0, 
ftFintra = 0.15 eV case. The diamonds are experimental data 
measured in K0.3M0O3 at T = 10 K 0. 



optical excitations. We determine the current-current 
correlation function, i.e. the optical conductivity of the 
related two-band model as a function of band filling. The 
transverse equation of motion approach has been used to 
collect the most singular contributions to the optical con- 
ductivity. It is shown that the present multiband optical 
analysis represents a generalization of the usual force- 
force correlation function method, and that in the DC 
limit it approaches correctly the results of Boltzmann 
equations, due to its gauge-invariant form. The present 
optical conductivity model gives the frequency and tem- 
perature dependences of the single-particle contributions 
to the optical conductivity spectra in the ordered CDW 
state which are consistent with both the experimental ob- 
servation and the prediction of the longitudinal response 
theory. It is explained that, while the response to the lon- 
gitudinal fields is associated with the direct electron-hole 
pair excitations, the response to the transverse electro- 
magnetic fields can be understood in terms of the indirect 
(impurity-assisted) electron-hole pair excitations. 



Acknowledgements 

This research was supported by the Croatian Ministry 
of Science and Technology under the project 0119-256. 
One of the authors (I.K.) would like to acknowledge the 
hospitality of the Institute of Physics of Complex Matter, 
Lausanne, where parts of this work were completed. 

APPENDIX A: CURRENT AND STATIC 
RAMAN VERTICES 

The vertex functions in the expression (|25[1 depend on 
the unperturbed vertices J„(k) = (e/h)dei(k)/dk a and 
7i' Q (k; 2) = {m/h 2 )d 2 ei{k)/dkl in the way 

J Q LL '(k) = £l/ k (^)^(Z,LVa(k), 
l 

7 ^'(k;2) = ^L/ k (/,L)^(/,L') 7 L(k;2). (Al) 



For the a = a polarization of the electromagnetic fields, 
the result is 

jF' CC (k) = TK(k)-, 2 (k)]jr(k) 

= T cos^(k)J a cc (k), (A2) 

j£ c (k) = 2u(k)«(k)JS°(k) 

= sin^(k)J a cc (k), (A2') 

7 £p 7C (k;2) = T K(k)- V 2 (k)] 7 -(k;2) 

= Tcos^(k) 7 -(k;2). (A3) 

Similarly, for a = b 



jp cc (k) = j h -(k), jp(k) = 0, 

7^' CC (*2) = 7 -(k;2). 



(A4) 



Finally, using Eqs. flTSJl. ljA2() and (|A4|) we can check the 
Ward identity [3 which connects the intraband current 
vertex J^ L (k) with the electron group velocity v„(k) 

J^(k) = «£(k) = |^). (A5) 
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